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ABSTRACT  OF  THESIS 


A  MODELING  STUDY  OF  THE  CLOUD-TOPPED  MARINE  BOUNDARY  LAYER 


coupled  convective-radiative,  boundary-layer  model  of  marine  stratocumulus  clouds 


is  presented.  The  model,  which  slightly  generalizes  Lilly’s  (1968)  cloud- topped  mixed-layer 
model,  has  as  dependent  variables  the  cloud-top  height,  the  cloud-base  height,  mixed- 
layer  equivalent  potential  temperature  and  total  water  mixing  ratio,  the  turbulent  fluxes 
of  equivalent  potential  temperature,  total  water  mixing  ratio,  and  virtual  potential  tem¬ 
perature,  the  cloud-top  jumps  of  equivalent  potential  temperature  and  total  water  mixing 
ratio,  the  cloud-top  temperature,  and  the  net  radiative  flux  divergence  at  cloud  top  and 
in  the  mixed  layer. 


The  model  has  been  programmed  for  both  the  steady-state  and  time  dependent  cases 
using  two  different  closure  assumptions.  Closure  method  one  is  a  weighted  average  of 
Lilly’s  (1968)  maximum  and  minimum  entrainment  assumptions.  Closure  method  two 
assumes  the  ratio  of  the  integral  buoyant  energy  dissipation  to  the  buoyant  energy  pro¬ 
duction  to  be  constant.  Experiments  were  performed  to  determine  which  closure  method 
produced  the  best  results,  and  to  simulate  the  observed  diurnal  thinning  of  the  marine 
\  stratocumulus  clouds. 

Results  indicate  that  closure  method  two  provides  a  more  realistic  representation  of 
the  observed  marine  boundary  layer  than  does  closure  method  one  for  both  the  steady- 
I  state  case  and  the  time  dependent  case.  Results  from  the  diurnal  cycle  simulations  indicate 
I  that  it  is  possible  to  simulate  the  observed  thinning  trend  of  the  stratocumulus  cloud  layer 

(although  the  magnitude  of  the  thinning  is  less  than  observed.  Results  also  indicate  that 
t  thin  clouds  are  more  susceptible  to  thinning  than  are  thick  clouds, 
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Thomas  A.  Guinn 

Department  of  Atmospheric  Science 

Colorado  State  University 

Fort  Collins,  Colorado  80523 
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Chapter  1 


INTRODUCTION 

In  a  field  which  gains  the  majority  of  its  notoriety  from  the  prediction  of  life  threat¬ 
ening  weather  such  as  hurricanes,  tornadoes,  and  lightning  it  seems  odd  that  marine 
stratocumulus  clouds  would  need  to  be  studied.  However,  when  one  considers  that  the 
globally  averaged  stratiform  cloud  coverage  is  on  the  order  of  12-15%  (Randall,  1985a), 
their  impact  on  the  planetary  albedo  and  therefore  the  global  climate  becomes  of  great 
importance.  It  has  even  been  calculated  that  a  4%  increase  in  the  global  coverage  of  the 
earth  by  stratocumulus  type  clouds  could  actually  offset  the  the  2-3  K  predicted  rise  in 
global  temperature  stemming  from  the  carbon  dioxide  greenhouse  effect  (Randall  et  a/., 
1984).  The  role  of  stratocumulus  clouds  in  the  climate  of  our  planet  does  not  end  with 
radiative  effects.  Stratocumulus  clouds  also  play  a  major  role  in  the  turbulent  heat  and 
moisture  exchange  between  the  boundary  layer  and  the  free  atmosphere.  Clearly,  then, 
if  the  global  climate  is  to  be  accurately  simulated  an  understanding  of  marine  stratiform 
clouds  must  be  acquired. 

Marine  stratocumulus  are  typically  found  over  the  cool  waters  of  the  subtropic  oceans 
where  large  scale  sinking  motion  suppresses  strong  convection,  near  the  eastern  edges  of 
the  subtropical  highs.  The  five  principal  areas  where  marine  stratocumulus  exist  are 
coastal  Oregon,  California  and  Mexico,  coastal  Equador,  Peru  and  Chile,  coastal  north¬ 
western  Africa,  coastal  southwestern  Africa,  and  coastal  northeastern  Africa  (Schubert  et 
af.,  1979). 

Because  of  their  widespread  global  coverage  and  their  impact  on  the  global  climate, 
it  is  no  wonder  that  extensive  modeling  studies  of  marine  stratocumulus  clouds  have  been 
undertaken,  the  goals  of  which  are  to  better  understand  the  interaction  of  the  variables 
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and  processes  which  govern  the  formation,  maintenance,  and  dissipation  of  these  clouds. 
The  difficulty  in  such  modeling  comes  from  the  extensive  parameterization  of  such  complex 
physical  processes  as  radiation,  entrainment,  and  turbulent  fluxes. 

The  first  to  develop  a  model  of  the  cloud-topped  marine  boundary  layer  was  Lilly 
(1968).  His  model  emphasized  longwave  cooling  off  cloud  top  as  the  important  driving 
mechanism  for  stratocumulus.  This  elegant,  one-dimensional,  mixed  layer  model  was 
soon  followed  by  several  models  which  attempted  to  generalize  his  mixed  layer  theory 
(Schubert,  1976;  Kraus  and  Schaller,  1978a;  Schubert  et  a !.,  1979).  More  complicated, 
higher-order  models  of  the  marine  stratocumulus  clouds  have  also  been  developed  which 
are  able  to  represent  the  marine  boundary  layer  with  significant  success  (e.g.,  Bougeault, 
1985).  However,  the  accuracy  of  the  simpler,  more  computationally  efficient  mixed-layer 
models  suggest  that  they  will  continue  to  play  an  important  role  in  understanding  the 
marine  boundary  layer.  Randall  (1985b)  best  described  the  importance  of  mixed  layer 
models  when  he  said,  “Although  complex  models  are  capable  of  simulating  many  details  of 
the  cloud-topped  marine  boundary  layer  with  remarkable  accuracy,  simple  models  continue 
to  be  of  great  importance,  for  reasons  that  go  beyond  their  computational  efficiency.  It 
is  often  more  challenging  to  represent  a  process  in  a  simple  model  than  in  a  complex  one; 
the  realistic  incorporation  of  a  process  into  a  simple  model  marks  the  achievement  of  real 
understanding.”  For  these  reasons,  we  have  chosen  to  study  the  marine  boundary  layer 
and  the  effect  of  the  diurnal  cycle  with  a  simple  one-dimensional  mixed  layer  model. 

The  model  presented  here  is  a  slight  generalization  of  Lilly’s  ( 1968)  mixed  layer  model. 
For  this  study,  Lilly’s  (1968)  model  is  generalized  by  considering  both  the  net  radiative 
flux  divergence  at  cloud  top  as  well  as  the  net  radiative  flux  divergence  within  the  mixed 
layer  to  be  important.  We  also  allow  the  absorbed  shortwave  radiation  to  vary  with  the 
sun  angle  to  enable  us  to  simulate  the  diurnal  cycle.  In  addition  to  these  changes,  Lilly’s 
(1968)  closure  equation  is  slightly  modified.  We  have  chosen  two  closure  assumptions 
for  this  study.  The  first  closure  is  a  weighted  average  of  Lilly’s  (1968)  maximum  and 
minimum  entrainment  case.  This  method  was  first  proposed  by  Schubert  (1976).  The 
second  closure  method  considers  the  ratio  of  the  integral  buoyant  energy  dissipation  to 
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the  integral  buoyant  energy  production  to  be  constant.  The  models  for  both  closure 
methods  are  presented  for  the  steady  state  case  as  well  as  the  time  dependent  case. 

The  outline  of  this  paper  is  as  follows.  Chapter  two  is  a  review  of  Lilly’s  original 
mixed-layer  model  while  chapter  three  describes  the  revisions  and  generalizations  made  to 
his  model  and  the  radiation  parameterization.  Experimental  results  of  our  models  for  both 
the  steady-state  and  time  dependent  cases  as  well  as  the  observational  studies  of  Betts 
(1988)  and  Davies  and  Blaskovic  (1988)  are  presented  in  chapter  four.  These  experiments 
focus  mainly  on  a  comparison  of  the  two  closure  assumptions  and  the  simulation  of  the 
observed  diurnal  thinning  of  stratocumulus  cloud  decks.  Finally,  our  conclusions  are 
presented  in  chapter  five. 


Chapter  2 


REVIEW  OF  LILLY’S  CLOUD  TOPPED  MIXED  LAYER  MODEL 

In  a  paper  which  has  become  a  classic  in  the  study  of  the  marine  boundary  layer, 
Lilly  (1968)  presents  two  boundary  layer  models — one  for  the  dry  boundary  layer  case  and 
one  for  the  cloud-topped  case.  This  chapter  is  devoted  to  reviewing  the  cloud-topped  case 
of  Lilly’s  paper. 

Before  looking  at  the  mechanics  of  Lilly’s  model,  it  is  first  necessary  to  list  the  as¬ 
sumptions  that  are  built  in.  The  first  assumption  is  that  shear  generated  turbulence  is 
confined  to  the  lowest  few  tens  of  meters  and  can  therefore  be  neglected.  The  radiation 
off  the  cloud  top,  however,  is  considered  to  be  essential  (Lilly,  1968)  and  strong  enough  to 
thoroughly  mix  the  boundary  layer  through  turbulent  processes.  Secondly,  the  air  above 
the  mixed  layer  (boundary  layer)  is  assumed  to  be  non-turbulent  and  of  high  enough 
temperature  that  if  forced  into  the  mixed  layer,  evaporative  cooling  would  not  lower  its 
temperature  beyond  that  of  the  mixed  layer.  If  this  happens,  the  air  at  cloud  top  would 
become  negatively  buoyant  and  penetrate  freely  into  the  mixed  layer  until  the  cloud  was 
completely  evaporated.  For  a  change  in  mixing  ratio  across  the  inversion  of  5-10  g/kg, 
an  increase  in  potential  temperature  of  ~  13  -  26°C  is  required  for  no  change  in  tempera¬ 
ture  to  occur  upon  evaporative  mixing  (Lilly,  1968).  In  more  concise  terms,  Lilly’s  (1968) 
condition  for  a  stable  cloud  layer  is  that  the  equivalent  potential  temperature  of  the  air 
above  the  inversion  must  be  at  least  as  great  as  that  of  the  air  in  the  mixed  layer.  Thirdly, 
the  mixed  layer  is  assumed  to  be  free  of  horizontal  variations,  i.e.,  it  is  assumed  to  be 
horizontally  homogeneous.  Finally,  the  large  scale  horizontal  divergence  and  surface  data 
are  assumed  to  be  known.  With  these  assumptions  in  hand,  Lilly’s  cloud-topped  model 
can  be  derived. 
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For  a  layer  in  which  no  condensation  takes  place,  the  potential  temperature  (9)  and 
water  vapor  mixing  ratio  (^)  can  be  regarded  as  conservative  thermodynamic  variables. 
In  addition,  if  the  layer  is  well  mixed,  9  and  q  can  be  regarded  constant  within  the  mixed 
layer.  In  a  cloud  topped  mixed  layer,  9  and  q  are  no  longer  appropriate  conservative 
thermodynamic  variables  because  condensation  must  now  occur.  However,  if  the  resulting 
liquid  water  (l)  is  not  precipitated  but  rather  is  carried  along  with  the  air,  the  total  water 
mixing  ratio  (q  +  l)  and  equivalent  potential  temperature  ( 9e )  can  then  be  considered  the 
appropriate  conservative  thermodynamic  variables.  Therefore,  in  a  well  mixed  boundary 
layer  of  depth  H,  9e  and  q  +  t  remain  constant  in  height.  It  should  be  noted  that  the 
omission  of  precipitation  restricts  the  model  to  rather  thin  cloud  layers. 

Assuming  the  mixed  layer  is  shallow  compared  to  the  scale  height  of  the  atmosphere, 
the  thermodynamics  can  be  simplified  through  the  use  of  a  linear  condensation  model. 
To  derive  this  relationship,  the  Clausius-Clapeyron  equation  is  combined  with  the  hydro¬ 
static  equation.  The  result  is  a  differential  relationship  between  saturation  mixing  ratio, 
potential  temperature,  and  height.  This  relationship  is  given  by 

dqaat  —  ad9  —  bdz,  (2.1) 

where 

Lq,at  (  P  \ 

°  RVT9  \p-ej  ’ 

and 

i  gpqnt  (  L  _  Cp  \ 

(p  _  et)  cp  \RvT2  RT ) 

(see  Appendix  for  the  complete  derivation  and  definitions  of  terms).  Although  a  and  b  are 
both  rapidly  changing  functions  of  height,  Lilly  maintains  that  for  a  shallow  layer  they 
may  be  considered  constant. 

The  shallow  layer  approximation  can  also  be  used  to  simplify  our  definition  of  9e. 
Lilly  (1968)  defines  9t  as  a  measure  of  the  total  sensible  plus  latent  heat  of  the  air.  For  a 
shallow  layer  where  T  «  9  ss  0e,  we  can  express  9e  as  a  differential  relationship  between 
potential  temperature  and  water  vapor  mixing  ratio.  Specifically, 

d9e  =  d9+  —dq. 

cp 


(2.2) 


With  the  use  of  (2.1)  and  (2.2),  one  can  easily  derive  an  expression  for  the  cloud 
base  height.  This  is  done  by  integrating  (2.1)  upwards  from  the  sea-surface  to  cloud  base 
( h ).  Potential  temperature  can  then  be  eliminated  from  the  resulting  expression  through 
substitution  of  (2.2).  The  end  result  is 


bh  =  (  1  +  a)[$s  -  (<?  +  01  ~  a(^s  -  *•) 


(2.3) 


where  qs  and  0es  are  the  saturation  mixing  ratio  and  saturation  equivalent  potential 
temperature  at  the  sea  surface,  and  a  =  (aL/cp)  where  L  is  latent  heat  and  cp  is  the 
specific  heat  of  dry  air  at  constant  pressure.  Here  and  throughout  this  paper,  the  subscript 
S  will  refer  to  those  properties  at  the  surface. 

Again  assuming  that  the  boundary  layer  depth  is  much  smaller  than  the  scale  height 
of  the  atmosphere,  it  can  be  shown  that  the  vertical  heat  flux  is  a  linear  function  of  height. 
This  is  done  by  first  combining  the  thermal  equation,  d9e/dt  =  0,  with  the  Boussinesq  form 
of  the  continuity  equation  to  obtain  the  thermal  equation  in  flux  form.  By  imposing  the 
assumption  of  horizontal  homogeneity  (i.e.,  dOe/dx  =  dOc/dy  =  0)  and  the  assumption 
that  the  boundary  layer  is  well  mixed  in  0t  (i.e.,  dOJdz  =  0)  the  thermal  equation  is 
reduced  to 


d9 e  d 

aT  ■  ~T,W- 


(2.4) 


Since  0e  is  constant  with  height,  it  follows  from  (2.4)  that  i u'0,e  is  a  linear  function  of 
height.  This  linear  relationship  can  be  expressed  as 


=  (l  -  §)  (WK)s  +  (2.5) 

where  the  subscript  H  refers  to  those  properties  just  below  cloud  top.  By  substituting 
(2.5)  into  (2.4),  the  thermal  equation  for  the  mixed  layer  is  written  in  final  form  as 

M.  =  [(Tgk  -  (g%)ffl_  (2  6) 

dt  H 


By  analogy,  starting  with  d(q  +  i)/dt  =  0,  the  total  water  mixing  ratio  equation  for  the 
mixed  layer  is  written  as 
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At  the  upper  boundary  (cloud  top),  the  thermal  equation  takes  a  different  form.  For 
cloud  top,  Lilly  (1968)  states  that  if  a  stable  temperature  discontinuity  is  to  exist,  it  must 
be  maintained  by  an  infinite  heating  or  cooling  rate  associated  with  a  discontinuity  in  the 
radiative  and/or  turbulent  heat  fluxes.  For  this  reason,  the  thermal  equation  at  cloud  top 
must  now  include  both  the  vertical  variation  in  de  and  the  radiative  flux  ( F ).  With  this 
in  mind,  the  thermal  equation  at  cloud  top  is  written  as 


dde 


d9e 


dF 


af— 


(2.8) 


If  (2.8)  is  integrated  from  just  below  cloud  base  ( H~ )  to  an  equal  distance  just  above 
cloud  base  (ff+)  the  following  expression  is  obtained 

rH+  doe 


J 

Jh- 


dt 


dz  =  -WH(6eU  ~  #e)  +  (u/0' )tf  ~  Fu 


(2.9) 


where  the  subscript  U  refers  to  those  properties  just  above  cloud  top.  However,  (2.9)  by 
can  be  simplified  by  first  using  Leibniz’  rule  to  expand  the  left  hand  side,  i.e., 


L 


d  rH+  dH 

=  -  /  Oedz-™(9eU-0e). 


H-  dt  dz  dt  Jh- 


dt 


(2.10) 


As  we  take  the  limit  A H  — *■  0,  the  first  term  on  the  right  hand  side  vanishes,  thus  the 
heat  balance  equation  at  cloud  top  can  be  expressed  as 

(^T  ~  A0e  +  (wf67)h  =  Fu,  (2.11) 

where  A0e  =  deu  -  0e  and  Fu  is  the  net  radiative  cooling  at  cloud  top.  By  analogy  it  can 
be  shown  that 

~  wh')  &{q  +  t)  +  +  i')H  =  0,  (2.12) 

where  A (q  +  t)  =  qu  ~(q  +  t). 

At  the  lower  boundary  (sea-surface),  w'9'e  can  be  related  to  the  near  surface  wind 
and  the  sea  surface  to  air  equivalent  potential  temperature  difference  through  the  use  of 
a  bulk  transfer  coefficient,  Cj-  By  analogy,  +  l')  is  related  to  the  near  surface  wind 
and  the  sea  surface  to  air  water  vapor  mixing  ratio  difference.  The  resulting  relationships 
are  given  as 


(«W.)s  =  CTV(9eS  -  9e) 


(2.13) 
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and 

^(T+lOs  =  CTV[qs  -  (q  +  l)}.  (2.14) 

The  wind  speed,  V,  is  typically  measured  from  ten  meters  above  the  surface. 

In  order  to  close  this  system  of  equations,  it  is  useful  to  discuss  the  concept  of  virtual 
potential  temperature  ( 9V ).  The  virtual  potential  temperature  is  defined  in  the  usual  way 
as  the  potential  temperature  a  parcel  of  air  would  have  if  it  had  the  same  density  as  dry 
air  at  the  same  pressure.  It  is  therefore  a  measure  of  the  relative  buoyancy  of  the  air. 
Lilly  (1968)  defines  9V  to  include  the  contributions  of  both  vapor  and  liquid  water  to  the 
air  density.  Specifically, 

d9v  =  d9  +  9(6dq  -  dt)  (2.15) 


where  6  =  [(1  —  e)/e]  and  €  is  the  ratio  of  the  molecular  weight  of  water  to  the  molecular 
of  dry  air.  The  value  of  6  used  here  is  0.608.  Since  (2.2)  and  (2.15)  are  valid  for  both 
mean  and  fluctuating  quantities,  they  can  be  expressed  in  terms  of  turbulent  fluxes,  i.e., 

=  ^  +  (2.16) 

cp 

=  riV  +  9(6^  -  ~EFV).  (2.17) 


Using  (2.16)  and  (2.17),4he  turbulent  virtual  potential  temperature  flux  can  be  ex¬ 
pressed  in  terms  of  more  conservative  quantities  by  eliminating  the  potential  temperature 
flux  between  (2.16)  and  (2.17).  The  resulting  expression  is 


i SOI  =  W9'e-9v/(q'  - Vt')- 


--9(6  +  1) 


ri/q'. 


(2.18) 


Although  (2.18)  adequately  describes  the  turbulent  virtual  potential  temperature  flux 
in  the  boundary  layer,  it  is  a  discontinuous  function  across  the  cloud  base  because  of 
the  liquid  water  term.  Beneath  the  cloud  base,  the  liquid  water  vanishes  allowing  some 
cancellation  between  the  third  and  fourth  terms  of  (2.18).  In  the  cloud,  the  mixing  ratio 
becomes  saturated,  therefore  (2.1)  evaluated  at  constant  height  can  be  combined  with 
(2.2)  to  show  that 

v/q>  =  -v/9g.  (2.19) 

(1  +  a) 
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Substituting  (2.19)  into  (2.18)  allows  the  second  and  fourth  terms  of  (2.18)  to  be  combined. 
Under  these  two  conditions,  (2.18)  can  be  rewritten  as 

w'9'e  -  (Jr  -  w'q',  for  0  <  z  <  h; 

(2.20) 

[1+ffio0  ~  *  +  *0.  for  h<z<H. 

Having  discussed  the  virtual  potential  temperature,  it  is  now  possible  to  return  to 
the  closure  equation.  Lilly  (1968)  proposes  from  considerations  of  the  turbulent  energy 
balance,  that  the  equation  necessary  to  close  this  system  of  equations  is  one  which  predicts 
the  rate  of  turbulent  entrainment.  He  offers  two  cases,  one  for  the  maximum  entrainment 
case  and  one  for  the  minimum  entrainment  case.  These  are 


Maximum  :  /  ( wl9lv)dz  =  0,  but  ^  0  somewhere. 

Jo 

_  _ 

Minimum  :  (u/0') min  =  0*  but  /  ( xu‘0'v)dz  >  0. 

Jo 


(2.21) 

(2.22) 


The  maximum  entrainment  case  implies  that  the  heat  transported  down  from  above 
the  inversion  equals  that  cafried  up  from  the  surface.  Viscous  dissipation  of  energy  is 
therefore  neglected.  The  minimum  case  confines  all  the  energy  dissipation  within  the  re¬ 
gion  of  positive  conversion  from  potential  energy.  No  energy  is  transported  to  do  work 
outside  the  boundaries  of  the  positive  conversion  region  (Lilly,  1968).  Since  viscous  dissi¬ 
pation  is  sometimes  large,  observations  show  that  the  boundary  layer  tends  more  towards 
the  minimum  entrainment  case  (Lilly,  1968). 


Chapter  3 


MODELING  THEORY 

Presented  in  this  chapter  are  both  the  time  dependent  and  steady-state,  horizontally 
homogeneous,  cloud- topped  mixed  layer  models  based  primarily  on  the  work  of  Lilly  (1968) 
and  Schubert  et  a 1.  (1979).  The  conservative  thermodynamic  variables  are  equivalent  po¬ 
tential  temperature  6e  and  total  water  mixing  ratio  q  +  i.  Some  of  the  differences  between 
these  and  Lilly’s  (1968)  model  are:  (1)  radiation  is  allowed  to  penetrate  into  the  boundary 
layer;  (2)  cloud  top  values  of  longwave  radiation,  equivalent  potential  temperature,  and 
water  vapor  mixing  ratio  are  linear  functions  of  height  derived  from  climatological  data  for 
California  coastal  stations;  (3)  two  closure  assumptions  are  explored,  the  first  assumes  a 
weighted  average  of  Lilly’s  (1968)  maximum  and  minimum  entrainment  theories,  the  sec¬ 
ond  assumes  the  ratio  of  the  integral  buoyant  energy  dissipation  over  the  integral  buoyant 
energy  production  is  constant. 

Section  3.1  describes  in  detail  the  revisions  made  to  Lilly’s  (1968)  model  and  the 
two  closure  assumptions.  Section  3.2  decribes  the  radiation  parameterization.  Section  3.3 
describes  the  combined  radiative-convective  model  for  both  closure  assumptions.  Finally, 
section  3.4  discusses  the  steady-state  solutions  to  the  combined  radiative-convective  model, 
again  for  both  closure  assumptions. 

3.1  REVISIONS  TO  LILLY’S  MODEL 

In  Lilly’s  (1968)  model,  the  ratdiation  is  confined  to  am  infinitesimally  thin  layer  at 
cloud  top.  Observationad  results  from  Slingo  et  a /.  (1982)  and  Albrecht  et  ai.  (1985), 
however,  indicate  that  substamtial  radiative  cooling  occurs  in  the  mixed  layer.  For  this 
reason,  the  models  presented  in  this  chapter  lift  this  restriction  of  longwave  cooling  to 
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cloud  top  and  thus  allow  radiation  to  penetrate  into  the  mixed  layer.  The  immediate 
result  of  this  assumption  is  the  inclusion  of  the  vertical  variation  of  the  radiative  flux  in 
the  mixed  layer  thermal  equation  (2.4).  The  new  mixed  layer  thermal  equation  is  written 
as 

f'  w 

where  F  is  the  radiative  flux.  Since  9e  is  well  mixed  in  the  vertical,  it  follows  from 
(3.1)  that  ( w'9'e  +  F)  is  a  linear  function  of  height.  Further,  if  the  radiative  flux  alone  is 
assumed  to  be  linear  with  height,  then  u/0/e  must  also  be  linear  with  height.  These  linear 
relationships  are  given  as 

F  =  (l  -  -|)  Fs  +  jjFH  (3.2) 

and  (2.5)  respectively.  By  substituting  (2.5)  and  (3.2)  into  (3.1),  the  mixed  layer  thermal 
equation  becomes 

dOe  _  (tt/AQs  ~  +  Fs  -  Fh 

dt  ~  H  ‘ 

The  cloud  top  thermal  equation  (2.11)  is  also  slightly  revised  The  new  cloud  top 
thermal  equation  takes  the  form 

{W  ~  WfI)  A*e  +  =  FV~F"  (3-4) 


As  can  be  seen,  (3.4)  is  the  same  as  (2.11)  with  the  exception  that  Fu  —  Fh  now  represents 
the  net  radiative  cooling  at  cloud  top.  This  distinction  is  made  for  consistency  with  (3.3). 

The  second  revision  made  to  Lilly’s  (1968)  model  is  in  the  calculation  of  the  cloud 
top  values  of  9e  and  q.  Instead  of  holding  these  variables  constant,  they  are  allowed  to 
vary  linearly  with  height.  The  linear  functions  used  in  this  model  were  derived  from  July 
averaged  sounding  data  for  San  Diego  and  Oakland,  CA  taken  from  the  U.S.  Dept,  of 
Commerce  Clim&tologic&l  Data,  National  Summary  for  the  years  1976-1980.  These  linear 
relationships  are 

f  5.50  -  0.00066  x  H,  for  San  Diego; 
qu  =  {  (3-5) 

l  5.54  -  0.00112  x  H ,  for  Oakland 

and 


(  313.95  +  0.00566  x  H ,  for  San  Diego; 

9eU  =  {  ,  (3-6) 

l  311.15  +  0.00341  x  H,  for  Oakland 
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where  the  units  of  qu  and  &u  are  g  kg-1  and  Kelvin  respectively. 

A  third  revision  to  Lilly’s  (1968)  model  is  the  parameterization  of  the  cloud  top 
vertical  velocity  in  terms  of  the  large  scale  horizontal  divergence.  Because  the  vertical 
velocity  at  cloud  top  is  difficult  to  measure,  it  is  preferred  to  parameterize  it  in  terms  of 
a  known  climatological  mean.  In  this  case,  the  large  scale  horizontal  divergence  is  used. 
The  parameterization  is  accomplished  by  integrating  the  Bousinesq  form  of  the  continuity 
equation  from  the  sea-surface  to  cloud  top.  Since  the  large  scale  horizontal  divergence  is 


constant  with  height  in  the  boundary  layer,  the  integration  is  simply 

( du  dv\  dw  , 

V<9z  +  dy)  Jo  2  ~  Jo  dz  2' 

Carrying  out  the  integration  gives 

(du  dv\ 

-  (.£  +  Ty )  H  = 

where  w(z  =  0)  =  0  since  the  ocean  surface  is  assumed  to  have  no  slope.  By  defining 
D  =  (du/dz  +  dv/dy ),  (3.4)  and  (2.12)  can  be  rewritten  as 


+  DH^j  A9e  +  Wh  =  Fu-F„  (3.7) 

and 

(“  +  DH )  A(q  + 1)  +  +  =  0,  (3.8) 

Equations  (3.7)  and  (3.8)  are  both  predictive  equations  for  H.  In  order  that  they 
predict  H  in  a  consistent  manner, 


A(q  +  l)(u/0'e)H  -  A 9e  u/(q'  +  1%  =  A (q  +  1){FV  -  F„).  (3.9) 


Equation  (3.9)  is  derived  by  multiplying  (3.7)  by  A(g-M)  and  (3.8)  by  -A0e  and  summing. 

The  revisions  to  Lilly’s  (1968)  model  are  now  complete  with  the  exception  of  the 
closure  equation.  Two  methods  of  deriving  a  closure  equation  are  discussed  in  the  following 
subsections.  The  first  is  taken  from  the  work  of  Schubert  (1976)  where  a  weighted  average 
of  Lilly’s  (1968)  maximum  and  minimum  entrainment  cases  is  considered.  The  second 
considers  the  ratio  of  the  integral  buoyant  energy  dissipation  to  the  integral  buoyant 
energy  production  to  be  constant,  as  suggested  by  Kraus  and  Schaller  (1978b). 
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3.1.1  CLOSURE  METHOD  ONE 


Schubert  (1976)  suggests  that  perhaps  a  better  means  of  examining  the  marine  bound¬ 
ary  layer  is  to  consider  a  weighted  average  of  Lilly’s  maximum  and  minimum  entrainment 
hypotheses.  Using  this  concept,  the  closure  equation  takes  the  form 


k 

H 


H  w'6'vdz  +  a-k)(w'0'v)mia 


=  0, 


(3.10) 


where  0  <  k  <  1  and  the  equation  for  the  turbulent  virtual  temperature  flux  is  given  by 
(2.20)  and  can  be  rewritten  as 

{w'&t  -  7u>'(g/  -(- 1'),  for  0  <  z  <  h\ 

(3.11) 

-  §w'(q'+t’),  for  h  <  z  <  H. 

where 

7  3  (£  - *)  “■* 0  3  |1+(i+°)131-  <3-i2) 

It  should  be  mentioned  that  turbulent  liquid  water  flux  is  zero  for  0  <  z  <  h;  however,  it 
is  included  in  (3.11)  to  simplify  future  derivations. 

As  mentioned  previously,  is  linear  with  height  throughout  the  mixed  layer  be¬ 
cause  0e  is  constant  with  height  in  the  mixed  layer.  The  same  holds  true  for  u/(g'  + 1') 
since  (q  +  l)  is  also  constant  with  height  in  the  mixed  layer.  This  linear  expression  is 
analogous  to  (2.5)  and  is  given  as 


W{<t  +  V)  =  (l  -  | +  l')s  +  jv/tf  +  t%  (3.13) 

Because  r&Ve  and  v/(q'  +  V)  are  linear  functions  of  height,  it  follows  from  (3.11)  that  u/0/„ 
must  also  be  a  linear  in  height  but  with  a  discontinuity  across  cloud  base.  One  might 
think  then  that  the  minimum  in  ti/0{,  could  occur  at  one  of  four  places,  at  the  surface,  just 
below  cloud  base,  just  above  cloud  base,  or  at  cloud  top.  However,  a  closer  examination 
of  (3.11)  shows  that  the  jump  in  v/9 'v  across  cloud  base  must  always  be  positive.  This 
eliminates  the  possibility  of  a  minimum  in  u/0(,  occurring  just  above  cloud  base  since 
the  value  just  below  cloud  base  will  always  be  less.  To  show  this,  use  (3.11)  to  take  the 
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difference  between  (tt/0{,)A+  and  where  h+  and  h~  refer  to  a  height  just  above 

and  just  below  cloud  base  respectively.  After  some  algebraic  manipulation  the  result  can 
be  written  as 

=  (»(6  +  1)  -  f )  -  w'(q'  +  t%)  (3.14) 

Because  the  air  at  cloud  base  is  saturated,  we  can  substitute  (2.19)  into  (3.14).  Making 
this  substitution  and  cancelling  terms  gives 

-  (yKh-  =  [f  -  W  + 1)]  &F)k+  (3-15) 

Since  (w‘l')h+  >  0,  (a/a)  >  0,  and  (a/a)  >  9(6  +  1),  the  left  hand  side  of  (3.15)  must 
be  positive.  The  only  remaining  possible  locations  for  the  minimum  turbulent  virtual 
potential  temperature  flux  are  at  the  surface,  just  below  cloud  base,  or  at  cloud  top. 

Using  the  above  constraints,  the  integrated  closure  equation  can  be  derived  by  substi¬ 
tuting  (2.5),  (3.11),  and  (3.13)  into  (3.10)  and  integrating.  The  resulting  closure  equation 
relates  H ,  h,  (v/j%)s,  «7(7TF)s,  and  +  /')»•  This  relationship  is  written 


as 


[±T1  f2*  -  w)  + 1  “'<«'  + *>»  -  -  *)  (|j)  +  *  »'(«'  +  «')*} 


+  k 


the  minimum  of 


_  -  8  vStf  +  1% 

(i  -  ft)  ws + -  7  (i  -  ft)  + ns  - 

(w'O'Js  -  7  u>'(g'  +  f')s 


7ftty,(9/ +  f')w  J  =  0 

(3.16) 


3.1.2  CLOSURE  METHOD  TWO 


For  closure  method  two,  the  ideas  of  Kraus  and  Schaller  (1978b)  are  followed  in  that 
a  closure  equation  which  is  based  on  the  ratio  of  the  integral  buoyant  energy  dissipation 
over  the  integral  buoyant  energy  production  is  considered,  i.e., 


fH _  [H _ 

k  /  [u/0,v  >  0 ]dz  +  /  [u /0’v  <  0 ]dz  = 
Jo  Jo 


0. 


(3.17) 
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As  with  closure  method  one  however,  the  closure  equation  still  varies  with  the  shape 
of  the  virtual  potential  temperature  flux  profile.  It  is  therefore  assumed  that  (w'9'v)h  >  0 
and  only  the  three  profiles  given  in  Fig.  3.1  are  considered. 

The  first  step  in  determining  which  virtual  potential  temperature  flux  profile  to  use 
is  to  check  the  sign  of  (w'9'v)q.  If  it  is  positive,  Fig.  3.1a  is  the  appropriate  profile  to  use 
and  the  closure  equation  takes  the  form 


k  {  [(«>,0(()a+  +  {w'^h  (H  -  h)  +  2c(u^^v)s}  +  (h  -  zc)(w'6'v)h-  =  0  (3.18) 


where  zc  is  the  crossover  height  at  which  (w'9'v)  —  0  between  the  surface  and  cloud  base. 
This  crossover  height  is  given  by 


£ £  _  (w'Oy)s 

*  («w„)s  -  (y¥v)h. 


(3.19) 


Equations  (3.9),  (3.18),  and  (3.19)  form  a  nonlinear  system  of  equations  in  the  three 
unknowns  zc ,  (w,9‘e)fj,  and  u/(q'  +  V)h-  A  solution  is  found  by  algebraically  reducing  the 
system  of  equations  to  one  quadratic  equation  in  the  unknown  zc  and  iteratively  finding 
its  zero  using  a  simple  secant  algorithm. 

If  (t u'Q'v)s  <  0,  the  closure  equation  can  take  on  two  different  forms  corresponding  to 
Fig.s  3.1b  and  3.1c.  To  determine  a  solution,  it  is  first  assumed  that  ( w'O (,)/,-  <  0.  Under 
this  restriction,  the  closure  takes  the  form 


k  [(u/0^  +  (u/O* ]  (H-h)+  [(u^)*  +  K*' )*- 


h  =  0 


(3.20) 


Equations  (3.9)  and  (3.20)  form  a  linear  system  of  two  equations  in  the  unknowns  (u/0£)//> 
and  v/tf  +  V)n  which  can  easily  be  solved  algebraically.  Once  the  solution  is  found,  (3.11) 
may  be  used  to  determine  if  it  fits  the  assumed  profile  shape.  If  it  does  not,  the  shape 
of  the  virtual  potential  temperature  profile  is  assumed  to  be  that  of  Fig.  3.1c.  In  this 
case,  the  closure  equation  takes  the  form 


k  {(h  -  Ze){W9>v)h-  +  [(«/<?(,)„*  +  (u WV)H\  (H  -  /»)}  +  (W9'v)szc  =  0  (3.21) 


Equations  (3.9),  (3.19)  and  (3.21)  can  then  be  used  to  solve  for  {y7¥e)n  and  it/iq1  +  V)H 
in  the  same  manner  as  for  the  case  where  {w'9'v)s  >  0. 
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Figure  3.1:  The  shapes  of  the  three  possible  profiles  for  the  virtual  potential  temperature 
flux. 


17 


3.2  RADIATION  PARAMETERIZATION 

Since  longwave  cooling  off  cloud  top  is  the  important  driving  mechanism  in  the  ma¬ 
rine  boundary  layer,  the  radiation  paramaterization  warrants  a  detailed  description.  The 
philosophy  that,  consistent  with  the  limitation  of  the  vertical  thermal  structure  to  two 
degrees  of  freedom,  the  vertical  resolution  of  the  radiative  cooling  should  also  be  limited 
to  two  degrees  of  freedom  is  followed.  This  means  radiation  can  appear  at  most  in  both 
the  mixed  layer  thermal  equation  (3.3)  and  the  cloud  top  thermal  equation  (3.4).  This 
is  slightly  different  from  Lilly’s  (1968)  model  where  radiative  cooling  appears  only  in  the 
cloud  top  thermal  equation  (2.11). 

Schubert  (1979)  gives  the  equations  for  the  net  radiative  flux  as 

Fu-F„  =  (pcp)-1{(l  -  m)(*4  -  4)  -  (1  -  p')S},  (3.22) 

and 

Fh  ~  Fs  =  (pcP)-1  ~  4)  “  (3-23) 

where  p  is  the  air  density,  a  is  the  Stefan-Boltzmann  constant,  Th  is  the  cloud  top  temper¬ 
ature  in  Kelvin,  4  *s  downward  longwave  radiative  flux  in  Wm~2  attributed  to  that 
portion  of  the  atmosphere  which  lies  above  the  mixed  layer,  S  is  the  absorbed  broadband 
shortwave  radiative  flux  in  Wm'J,  and  p  and  p!  are  the  longwave  and  shortwave  radiation 
partitions  respectively.  These  partitions  can  take  on  values  from  zero  to  one.  If  both  are 
set  to  zero,  the  radiation  parameterization  reduces  to  Lilly’s  case. 

For  this  model,  the  downward  longwave  radiative  flux  term  in  (3.22)  and  (3.23)  is 
expressed  as  a  linear  function  of  cloud  top  height.  The  function  was  derived  based  on  July 
averaged  sounding  data  for  the  five  years  1976-1980  as  extracted  from  the  U.S.  Dept,  of 
Commerce  Climatological  Data,  National  Summary.  Two  California  locations  were  used, 
San  Diego  and  Oakland.  The  sounding  data  combined  with  the  midlatitude  ozone  profile 
as  taken  from  the  U.S.  Standard  Atmosphere  Table,  1976  and  a  uniform  carbon  dioxide 
profile  of  .501gkg_1  was  entered  into  a  broadband  longwave  radiation  model  described 
by  Cox  (1973).  The  model  output  was  used  to  calculate  the  linear  relationships  for  the 
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two  locations  which  are  strictly  functions  of  cloud  top  height.  The  resulting  equations  are 
given  as 

(  333.2  -  0.02980  x  H,  for  San  Diego; 

4  =  (3-24) 

l  333.2  -  0.02799  x  H,  for  Oakland 

where  C\j  is  in  units  of  Wm'2. 

The  calculation  for  S  was  made  by  entering  the  above  sounding  data  into  a  broadband 
shortwave  radiation  model  described  by  Manabe  and  Wetherald  (1967).  Since  S  varied 
little  with  height  at  the  levels  of  typical  marine  stratocumulus  cloud  tops,  it  is  considered 
to  be  constant  in  height.  The  S  calculation  is,  however,  allowed  to  vary  with  time.  This 
allows  for  the  inclusion  of  a  diurnal  cycle  into  the  model.  For  simplicity,  only  solar  orbital 
geometry  is  considered  in  calculating  the  diurnal  cycle.  In  short,  the  S  calculation  is 
taken  to  be  only  a  function  of  the  cosine  of  the  solar  zenith  angle.  The  dependence  of 
upon  S  optical  path  length  is  neglected.  Figure  3.2  shows  the  difference  between  S  as 
calculated  without  consideration  cl  the  optical  path  dependence  and  S  as  calculated  with 
consideration  the  optical  path  dependence.  As  can  be  seen,  the  accuracy  gained  is  not 
warranted  for  the  amount  of  complexity  that  must  be  added  to  include  the  optical  path 
dependence.  In  addition,  its  inclusion  would  require  addition*.  degrees  of  freedom  for  the 
radiation  parameterization,  which  as  mentioned  previously,  is  not  desirable. 

The  constraint  on  S  is  that  it  reaches  its  climatological  maximum  at  local  solar  noon 
and  goes  to  zero  at  both  sunrise  and  sunset.  The  maximum  value  for  S  (corresponding 
to  local  solar  noon)  used  for  this  model  was  taken  from  the  900mb  level  of  the  broadband 
shortwave  radiation  model  output.  This  value  was  then  multiplied  by  a  cloud  shortwave 
absorption  coefficient  of  7%.  This  absorption  coefficient  corresponds  to  Lilly’s  (1968)  value 
which  was  calculated  from  the  observational  results  of  Neiburger  (1949).  The  calculated 
maximum  values  of  S  are 

{75.2,  for  San  Diego; 

(3.25) 

73.5,  for  Oakland, 


where  Sm^x  is  in  units  of  Wm*1. 


SHORTWAVE  ABSORPTION  (Wm~2) 
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HOURS  (LST) 


Figure  3.2:  The  absorbed  broadband  shortwave  radiation.  The  dashed  line  represents 
calculations  using  solar  geometry  alone.  The  solid  line  represents  calculations  using  both 
solar  geometry  and  optical  path  dependence. 


20 


In  summary,  the  equation  for  S  is  written  as 

'  2 


=  *(;)  {si 


in  <f>  sin  d  +  cos  <j>  cos  d  cos  r} 


(3.26) 


where  <f>  is  latitude,  So  is  <5max  normalized  by  the  local  noon  value  of  the  cosine  of  solar 
zenith  angle  and  the  square  of  the  distance  function  (r/r),  F  and  r  being,  the  mean  and 
instantaneous  distance  of  the  earth  from  the  sun  respectively,  d  is  the  solar  declination 
angle,  and  r  is  the  solar  hour  angle,  i.e.  the  angular  distance  of  a  point  on  the  earth  from 
local  solar  noon.  Note  that  the  bracketed  term  in  (3.26)  is  simply  the  definition  of  the 
cosine  of  the  solar  zenith  angle.  July  averages  for  the  calculations  of  r  and  d  are  based  on 
data  extracted  from  the  Smithsonian  Meteorological  Tables,  6th  edition. 

Since  S,  r,  r,  <f>,  and  d  are  constant  with  time  (within  the  time  frame  of  a  few  days), 
all  time  dependence  in  the  shortwave  radiation  calculation  must  stem  from  the  hour  angle. 
This  time  dependence  can  be  written  as 


r  =  jj(«  -  12),  (3.27) 

where  r  is  in  radians  and  t  is  the  local  sun  time  in  hours. 

In  order  to  calculate  the  radiative  terms  ( Fu  -  Ffj)  and  (Fu  -  Fs ),  an  equation 
for  Th  is  needed.  This  equation  is  derived  first  starting  with  the  equation  for  potential 
temperature  at  cloud  top  (Off).  The  potential  temperature  at  cloud  top  is  equal  to  the 
potential  temperature  at  cloud  base  (0^)  plus  the  change  that  occurs  when  following  a 
moist  adiabat  from  cloud  base  to  cloud  top.  That  is 

e„  =  eK+(^)t(n-h).  (3.28) 

However,  (d0/dz)et  is  constant.  To  show  this,  simply  set  d9e  =  0  in  (2.2)  and  substitute 

into  (2.1).  The  result  is 

=  (329) 

In  addition,  0*  can  be  calculated  from  our  conservative  thermodynamic  variables  by  use 
of  (2.2),  i.e., 

9h  =  0e-  —(q  + 1). 
cp 


(3.30) 
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Substituting  (3.29)  and  (3.30)  into  (3.28)  while  making  use  of  Poisson’s  equation  gives 


Th  = 


9e  -  -( q  +  t)  +  Ci(H  -  h) 


(S' 


(3.31) 


where  C\  =  (6/a)[a/(l+a)],  po  is  1000  mb,  pjj  is  the  pressure  at  cloud  top,  and  k  =  ( R/cp ) 
with  R  and  cp  being,  respectively,  the  ideal  gas  constant  for  dry  air  and  the  specific  heat 
of  dry  air  at  constant  pressure. 

Equation  (3.31)  is  not  complete,  however,  since  a  way  to  calculate  ( ph/Po)k  is  still 
needed.  To  do  this,  it  is  convenient  to  express  the  hydrostatic  equation  in  terms  potential 
temperature.  The  resulting  expression  is 


9  \PoJ 


(3.32) 


where  g  is  gravity.  By  integrating  (3.32)  from  the  surface  to  cloud  base,  an  expression  for 
( Ph/Po)k  can  be  derived.  This  integration,  however,  must  be  done  in  two  steps  since  9  is 
constant  below  cloud  base  but  is  height  dependent  above  cloud  base. 

The  first  step  is  to  integrate  from  the  surface  to  cloud  base  and  the  second  step  is 
to  integrate  from  cloud  base  to  cloud  top.  Integration  of  (3.32)  from  the  surface  to  cloud 
base  gives 

1(9' -©']-£■ 

The  integration  of  (3.32)  from  cloud  base  to  cloud  top,  however,  is  slightly  more  difficult 
because  of  the  height  dependence  of  9.  The  equation  for  9{z)  to  use  in  this  case  is  a  more 
general  form  of  (3.28).  That  is 


9(z)  =  9h+Cl(z-h). 

Substituting  (3.34)  into  (3.32)  and  integrating  from  cloud  base  to  cloud  top  gives 

CijH  -  h) 


(3.34) 


1  + 


!)- 


[©-©>*'■(-  .. 

By  summing  (3.33)  with  (3.35)  we  obtain  the  desired  expression  for  (ptf/po)*,  i-e.. 


(3.35) 
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By  substituting  (3.36)  into  (3.31),  Th  can  be  expressed  in  terms  of  predicted  model 
variables  and  known  constants.  This  expression  takes  the  form 


T,,-[»*+C,(ff-4)]{(£)“-^--£rta(l  + 

(\po/  cp0/,  CpGi  v 


Ci(H-h)' 
Ok  . 


(3.37) 

In  summary,  with  the  use  of  (3.24),  (3.26),  and  (3.37),  the  radiation  variables  Fu  -  FH 
and  Ffj  —  Fs  can  be  expressed  completely  in  terms  of  predicted  model  variables  and  known 
constants. 


3.3  THE  COMBINED  CONVECTIVE-RADIATIVE  MODEL 


With  the  revisions  to  Lilly’s  (1968)  model  complete  as  well  as  the  radiation  param¬ 
eterization,  a  combined  convective-radiative  model  for  both  closure  methods  can  now  be 
written.  These  two  models  are  presented  in  the  following  two  subsections. 

3.3.1  CLOSURE  METHOD  ONE 

Closure  method  one  again  corresponds  to  the  weighted  average  of  Lilly’s  maximum 
and  minimum  entrainment  assumption.  The  model  equations  consist  of  the  surface  flux 
equations  (2.13)  and  (2.14),  the  cloud  base  equation  (2.3),  the  cloud  top  jump  definitions 
as  defined  in  (2.11)  and  (2.12),  the  cloud  top  temperature  equation  (3.37),  the  radiation 
equations  (3.22)  and  (3.23), "the  consistency  relation  (3.9),  the  mixed  layer  thermal  and 
total  water  budget  equations  (3.3)  and  (2.7)  respectively,  the  cloud  top  thermal  equa¬ 
tion  (3.4),  and  the  closure  equation  (3.16).  These  thirteen  equations  form  a  closed  set 
in  thirteen  unknowns.  Since  the  model  contains  three  predictive  variables  (0e,  (q  +  i), 
and  H ),  three  initial  conditions  are  needed  to  start  the  model.  It  should  be  noted  that 
any  reasonable  initial  conditions  of  the  predictive  variables  may  be  used  as  long  as  the 
model  is  run  long  enough  for  all  variables  to  begin  to  repeat  themselves.  For  the  sake 
of  computational  efficiency,  however,  the  initial  conditions  are  obtained  from  the  steady 
state  solutions  (which  will  be  discussed  in  the  following  section).  This  assures  the  user  of 
a  rapid  convergence  to  the  real  solution. 

With  the  initial  conditions  specified,  the  numerical  time  integration  is  performed  in 


the  following  order: 
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(w>d'e)s  =  CTV(0eS  -  Be) 

(3.38) 

t</( q '  +  F)s  =  CrV[qs  -  (q  +  £)] 

(3.39) 

l  (1  +  a)[?5  -  (q  +  *)]  -  a(Bes  ~  #e) 
h-  b 

(3.40) 

A0e  =  deU  -  de 

(3.41) 

A (q  +  £)  =  qu  -  (?  +  t) 

(3.42) 

Th  =  +  CX{H  -  /»)] 


JL. 

cp&h 


(3.43) 


Fu-Fh  =  (pcp)-'{(  1  -  MoTb  -  4)  -  (1  -  p')S },  (3.44) 


FH-Fs  =  {pcp)~l{p{aTA„  -  4)  -  p'S}t  (3.45) 


Oil  ®12 

v 

a2i  ajj 

,&2. 

dBt  _  (v/O'Js  -  (w'flQtf  4-  fg  -  Fh 
dt  ~  H 


(3.46,  3.47) 


(3.48) 


d(q  +  t)  _  (w'q')s  -  u/W  +  tf) h 

dt  ~  H 


(3.49) 


dH_  _  Fb-FH-WH  _  DH 
dt  A0e 
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With  the  exception  of  (3.46)  and  (3.47),  all  the  above  equations  are  taken  directly  from 
previous  text.  These  two  equations  are  simply  a  shorthand  matrix  notation  for  (3.9) 
and  (3.16).  The  dependent  variables,  externally  specified  parameters  and  constants  for 
the  above  equations  are  given  in  Tables  3. 1-3.4.  For  those  values  which  depend  on  a 
reference  temperature  and/or  pressure,  we  use  sea-surface  temperature  minus  3.5  degrees 
and  surface  pressure  minus  35  mb.  This  represents  an  approximate  midpoint  to  the 
boundary  layer. 

The  numerical  time  integration  is  best  demonstrated  by  the  following  six  step  iteration 
process.  These  six  steps  constitute  a  single  iteration  cycle  and  are  as  follows: 

1)  Use  the  initial  conditions  of  the  predictive  variables  to  calculate  the  surface  fluxes 
(vu'0'e)s  and  u/( q '  +  £')s  using  (3.38)  and  (3.39). 

2)  Calculate  the  cloud  base  height,  h,  using  (3.40). 

3)  Calculate  the  cloud-top  jumps  in  equivalent  potential  temperature  and  total  water 
mixing  ratio  using  (3.41)  and  (3.42). 

4)  Calculate  the  cloud-top  temperature  using  (3.43).  Using  this  value  of  Th,  calculate 
the  radiative  flux  divergences  using  (3.44)  and  (3.45). 

5)  Using  the  two-by-two  matrix  given  by  (3.46)  and  (3.47),  calculate  cloud  top  fluxes 
(v/PJh  and  u/tf  +  i')H. 

6)  Integrate  (3.48)-(3.50)  using  a  classical  fourth  order  Runge-Kutta  scheme  to  predict 
new  values  of  0e,  (q  +  i),  and  H.  When  completed,  return  to  step  one. 

The  above  six  step  iteration  is  all  straightforward  with  the  exception  of  step  five. 
Because  of  the  form  of  (3.16),  the  minimum  in  must  be  known  before  the  proper 
coefficient  matrix  and  the  column  vector  on  the  right  hand  side  of  (3.46)  and  (3.47)  can 
be  obtained.  This,  however,  requires  that  variables  we  are  trying  to  calculate,  ( w'Be)fj  and 
ti/( q '  +  be  known.  To  avoid  this  circular  reasoning,  the  minimum  in  v/0'v  is  assumed 

to  occur  just  below  cloud  base  and  both  {u/d'e)n  and  u/(g'  +  f')w  are  t^en  calculated. 
Equation  (3.11)  can  then  be  used  to  determine  if  the  minimum  in  actually  did  occur 
just  below  cloud  base.  If  so,  the  correct  solution  has  been  found.  If  not,  the  minimum  is 
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Table  3.1:  Dependent  Variables 


H 

h 

9. 

<1  +  1 
WK)s 

WK)h 

«W+T)S 

«w+nH 

WK)s 

W%)k- 

WK)h+ 

W)|f 

A9, 

A(?  +  0 
Th 

Fu  -  Fh 
Fh  ~  Fs 


cloud-top  height 
cloud-base  height 

mixed  layer  equivalent  potential  temperature 

mixed  layer  total  water  mixing  ratio 

surface  equivalent  potential  temperature  flux 

cloud- top  equivalent  potential  temperature  flux 

surface  total  water  flux 

cloud-top  total  water  flux 

virtual  potential  temperature  flux  at  the  surface 

virtual  potential  temperature  flux  just  below  cloud  base 

virtual  potential  temperature  flux  just  above  cloud  base 

virtual  potential  temperature  flux  at  cloud  top 

cloud  top  equivalent  potential  temperature  jump 

cloud  top  total  water  jump 

cloud  top  temperature 

cloud  top  radiative  flux  divergence 

mixed  layer  radiative  flux  divergence 
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Table  3.2:  Externally  Specified  Parameters 

D  large-scale  divergence 
V  wind  speed 

&eS  equivalent  potential  temperature  at  sea  surface  temperature  and  pressure 
qs  saturation  mixing  ratio  at  sea  surface  temperature  and  pressure 
9eu  equivalent  potential  temperature  just  above  cloud  top 
qt,  total  water  mixing  ratio  just  above  cloud  top 

C  U  downward  longwave  radiation  flux  just  above  cloud  top 
S  total  shortwave  radiation  absorption 
H  longwave  radiatibn  partition 
m'  shortwave  radiation  partition 
k  entrainment  parameter 


cp  specific  heat  at  constant  pressure 
Ct  bulk  transfer  coefficient 
g  gravity 

a  Stefan- Boltzman  constant 

6  see  (2.15) 

Table  3.4:  Constants  which  depend  on  reference  temperature  and/or  pressure 

p  density 
0  see  (3.12) 

7  see  (3.12) 

a  see  (2.1) 

b  see  (2.1) 

§  reference  potential  temperature 
C\  see  (3.31) 


L  latent  heat  of  condensation 
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assumed  to  occur  at  the  surface  and  the  above  procedure  is  repeated.  If  no  true  minimum 
has  yet  been  found,  the  procedure  is  repeated  for  a  minimum  at  cloud  top.  If  one  and  only 
one  solution  leads  to  no  contradiction,  the  true  solution  is  assumed  to  have  been  found. 

This  model  has  been  programmed  in  FORTRAN  and  will  run  interactively  on  an  IBM- 
compatible  personal  computer.  This  particular  version  of  the  program  is  titled  STRATUS 
Version  2.0.  STRATUS  2.0  allows  the  user  to  specify  the  sea-surface  temperature,  the 
large  scale  horizontal  divergence,  the  wind  speed,  the  entrainment  parameter,  the  short¬ 
wave  partition,  the  long  wave  partition,  and  the  geographical  location  (San  Diego,  CA  or 
Oakland,  CA).  Output  for  the  tenth  model  day  is  then  written  to  two  data  files,  STRA¬ 
TUS. DAT  and  FLUXES.DAT.  STRATUS.DAT  contains  a  time  series  of  the  mixed  layer 
equivalent  potential  temperature,  the  mixed  layer  total  water  mixing  ratio,  cloud  base, 
cloud  top,  and  a  fL6  indicating  the  location  of  the  minimum  of  the  virtual  potential 
temperature  flux,  r  LUXES.DAT  contains  a  time  series  of  the  equivalent  potential  tem¬ 
perature  flux  and  the  total  water  mixing  ratio  flux  at  both  cloud  top  and  the  surface,  and 
the  radiative  flux  divergence  at  both  cloud  top  and  in  the  mixed  layer. 

3.3.2  CLOSURE  METHOD  TWO 

As  before,  closure  method  two  assumes  the  ratio  of  the  integral  buoyant  energy  dissi¬ 
pation  to  the  integral  buoyant  energy  production  is  a  constant.  This  difference  affects  only 
equations  (3.46)  and  (3.47).  All  other  equations  remain  the  same  as  in  closure  method 
one.  Similarly,  the  numerical  time  integration  is  also  identical  to  that  of  closure  method 
one  with  the  exception  that  step  five  now  follows  the  iterative  process  discussed  in  section 
3.1.2. 

This  model  has  also  been  programmed  in  FORTRAN  and  will  run  interactively  on 
any  IBM-compatible  personal  computer.  The  program  is  again  titled  STRATUS  but  with 
a  version  number  of  2.1.  The  input  and  output  of  versions  2.0  and  2.1  are  the  same 
with  the  exception  of  the  flag  in  STRATUS.DAT.  This  flag  no  longer  indicates  where 
the  minimum  in  the  virtual  potential  temperature  flux  is  located,  instead  it  indicates  the 
shape  of  virtual  potential  temperature  flux  profile.  It  should  be  noted  that  if  the  solutions 
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for  (Vtf'e)#  and  w'(q'  -f  t')H  do  not  agree  with  any  of  the  three  profile  shapes  given  in 
Fig.  3.1,  the  program  will  print  an  error  message  and  abort.  So  far  this  problem  has  not 
been  encountered  for  physically  reasonable  input  data. 

3.4  THE  STEADY-STATE  SOLUTIONS 

Solving  the  model  for  the  steady-state  case  is  obviously  much  different  than  for  the 
time  dependent  case.  For  the  steady-state  case,  all  derivatives  with  respect  to  time  are  set 
to  zero  and  a  system  of  nonlinear  albebraic  equations  is  obtained.  The  method  used  to 
solve  for  the  steady-state  case  is  to  reduce  the  system  of  equations  down  to  one  equation 
in  H  and  use  a  simple  secant  method  algorithm  to  find  its  zero.  Presented  here  are  the 
steady  state  solutions  for  both  closure  methods. 

before  describing  the  steady-state  model  equations,  two  items  must  be  discussed. 
First,  the  diurnal  cycle  can  no  longer  be  defined  for  the  steady-state  case.  Instead  a  daily 
averaged  solar  radiative  flux  (S)  value  of  22.3  Wm-2  as  suggested  by  Lilly  (1968)  will  be 
used.  This  value  corresponds  to  approximately  1.5  hours  after  sunrise  or  1.5  hours  before 
sunset  (0615  hrs  and  1745  hrs  local  sun  time  respectively).  Second,  from  (3.49)  we  can  see 
directly  that  under  steady-state  conditions,  the  total  water  flux,  tv,(q/  +  l'),  is  constant 
with  height.  For  this  reason, 'all  subscripts  on  this  variable  are  dropped.  With  these  two 
thoughts  in  mind,  the  derivation  of  the  steady  state  solutions  can  now  be  presented. 

The  steady-state  solutions  are  begun  by  first  deriving  expressions  for  the  surface  flux 
of  equivalent  potential  temperature  and  total  water  flux  in  terms  of  constants  and  exter¬ 
nally  specified  parameters  (see  Tables  3. 2-3.4).  The  former  is  accomplished  by  eliminating 
the  dependent  variables  9e  and  {w,9,e)n  between  the  steady-state  forms  of  (3.38),  (3.48) 
and  (3.50).  The  resulting  expression  is 

(F<J  -  Fh)  -  (fs  -  Fh )  -  DH(0,v  -  e.s ) 

(  ,)s  ■ - 


(3.51) 
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The  latter  expression  can  be  derived  by  eliminating  ( q  +  t)  and  v/(qf  4-  P)H  between 
equations  (3.8),  (3.39),  and  (3.49).  This  expression  is  given  as 


w'itf  +  (■')  = 


DH(qs  -  gy) 

(»+«W  ' 


(3.52) 


The  mixed  layer  equivalent  potential  temperature  can  also  be  eliminated  from  our  cloud 
base  equation.  This  is  accomplished  by  substituting  (3.38)  and  (3.39)  into  (3.40).  The 
resulting  expression  is 


_  (1  +  a)w'(<f  -f  £')  -  q(u/fl')s 
bCTV 


(3.53) 


To  derive  an  expression  for  (u/ #')//,  9e  is  eliminated  from  (3.50)  by  substituting  in  (3.38). 
This  equation  is 


(w'Wff  =  (Fu  ~  F„)  -  — (t^Js  ~  DH(9eU  -  9eS )  (3.54) 

All  the  necessary  information  is  now  available  to  calculate  the  virtual  potential  tem¬ 
perature  flux  at  the  surface,  just  below  cloud  base,  just  above  cloud  base,  and  at  cloud 
top.  These  expressions  are  taken  directly  from  (3.11)  and  are  given  respectively  as 


(y&i)s  =  (w'^)s  -  + f), 

(3.55) 

WK)h-  =  (w'®'e)fc  -  ru’W  +  <')> 

(3.56) 

(&K)k*  =  -  9  w,(q/  +  V), 

(3.57) 

{*Vv)H  ~  -9v/(q’  +  P), 

(3.58) 

where  (it/0g)*  is  defined  using  (2.5)  as 

(yZh  =  (i  - (^)s  +  j;(yve)H.  (3.59) 


The  importance  of  equations  (3.55)-(3.59)  will  become  apparent  in  the  following  two 
subsections. 
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3.4.1  CLOSURE  METHOD  ONE 

As  before,  closure  method  one  corresponds  to  a  weighted  average  of  Lilly’s  (1968) 
maximum  and  minimum  entrainment  case.  The  closure  equation  for  this  case  is  derived 
by  simply  integrating  (3.10).  This  integration  is  easily  accomplished  since  u/,0'v  is  linear 
with  height.  The  difficulty  stems  from  not  knowing  the  location  of  Therefore, 

to  solve  this  problem,  it  is  assumed  that  (itf'^min  is  located  just  below  cloud  base.  Under 
this  restriction,  the  closure  equation  then  takes  the  form 

2-1-fc  —  {VO'Jk-  +  (w'0'v)H  +  (w'9'v)h+  +  jj[(w'd'v)s  +  (w'd'v )k-  -  (t u'9'v)h  -  (v70i)h+  ]  =  0. 

(3.60) 

Since  all  the  variables  except  H  are  known,  (3.60)  can  be  solved  using  a  simple  secant 
method  algorithm.  Once  the  solution  is  found,  it  must  be  checked  to  determine  if  the 
minimum  did  actually  occur  just  below  cloud  base.  If  it  does,  the  solution  is  accepted.  If 
it  does  not,  (u/0(,)min  is  assumed  to  be  located  at  the  surface  and  H  is  solved  for  using 
the  following  closure  equation. 

k\wr9[,)s  +  {wWJh  +  (wrK)h+  +  Jj[(vrK)s  +  ~  (w%)h  ~  (uTBDa+J  =  o. 

(3.61) 

Again  the  solution  is  checked  to  determine  if  the  closure  assumption  is  correct.  If  it  is, 
the  solution  is  accepted.  If  not,  it  is  assumed  that  no  solution  exists.  It  should  be  noted 
that  the  possibility  of  (u/0J,)min  occurring  at  cloud  top  has  been  neglected.  It  is  seen  from 
(3.54)  and  (3.58)  that  for  this  case  to  occur,  the  radiative  cooling  at  cloud  top  would  have 

9 

to  be  negligible.  Since  radiative  cooling  is  the  primary  driving  mechanism  for  this  model, 
it  is  unlikely  that  this  case  could  physically  occur. 

With  the  closure  equation  determined,  a  closed  set  of  equations  (3.43)-(3.45)  and 
(3.51)-(3.61)  in  the  unknowns  H,  h,  Th,  ( Fu-Fh ),  (Fh~Fs),  (t (w'9'e)h,  {w'K)s, 
w/(9#  + 1'),  (w'#')^,  (u >'0'v)h-,  and  (u/0(,)s  is  obtained.  The  solution  process 

can  best  be  explained  by  the  following  four  step  iteration  sequence. 

1)  Make  an  initial  guess  of  the  cloud  top  height  H  and  the  radiative  flux  divergences 


( Fu  -  Fff)  and  ( Fh  -  Fs ). 
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2)  With  the  above  initial  guesses,  use  (3.51)-(3.59)  to  calculate  in  order  ( w'0'e)s , 

w'(q'  +  i'),  h ,  (u>'0£)j/,  (v/d'v)s,  {w'6'v)h+,  and  ( vj'9'v)h .  In  addition,  calculate 

the  radiation  variables  Th ,  (Fu  -  Fu),  and  (Ffj  -  Fs )  using  (3.43)-(3.45). 

3)  Using  the  above  information,  check  to  see  if  (3.60)  is  satisfied  within  some  tolerable 
limit.  If  it  is  not,  use  the  secant  method  to  predict  a  new  H  and  return  to  step  one  until 

(3.60)  is  satisfied. 

4)  Once  a  solution  has  been  found,  check  to  see  if  (w'&Dmin  occurs  at  the  assumed  lo¬ 
cation.  If  it  does,  the  solution  has  been  found.  If  it  does  not,  repeat  steps  1-3  substituting 

(3.61)  in  for  (3.60). 

This  model  has  been  programmed  in  FORTRAN  and  will  run  interactively  on  any 
IBM-compatible  personal  computer.  The  program,  STRATUS  Version  1.0,  requires  the 
following  input:  sea-surface  temperature,  large  scale  horizontal  divergence,  surface  wind 
speed,  longwave  and  shortwave  radiation  partitions,  the  entrainment  parameter,  the  ge¬ 
ographic  location,  and  the  initial  guess  for  cloud  top.  Output  consists  of  a  table  printed 
to  screen  containing  values  of  the  mixed  layer  equivalent  potential  temperature,  mixed 
layer  total  water  mixing  ratio,  the  convective  total  water  flux,  the  equivalent  potential 
temperature  flux  at  cloud  top  and  the  surface  ,  the  virtual  potential  temperature  flux  at 
cloud  top,  just  above  cloud  base,  just  below  cloud  base,  and  at  the  surface,  the  radiative 
flux  divergence  at  cloud  top  and  in  the  mixed  layer,  the  temperature  at  cloud  top,  the 
cloud  top  height,  and  the  cloud  base  height. 

3.4.2  CLOSURE  METHOD  TWO 

As  before,  closure  method  two  is  based  on  the  idea  that  the  ratio  of  the  integral  buoy¬ 
ant  energy  dissipation  over  the  integral  energy  production  is  constant.  As  with  closure 
method  one,  the  possibility  of  multiple  closure  equations  based  on  different  virtual  poten¬ 
tial  temperature  flux  profiles  still  exists.  Therefore,  a  similar  procedure  to  that  of  closure 
method  one  is  used.  The  profile  shape  is  assumed  to  be  that  of  Fig.  3.1a  and  (3.18)  is 
used  as  our  closure  equation.  Once  the  solution  has  been  found,  it  is  checked  to  see  if  it 
satisfies  the  assumed  profile  shape.  If  they  match,  the  true  solution  is  assumed  to  have 
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been  found.  If  not,  the  profile  shape  is  assumed  to  be  that  of  of  Fig.  3.1b  and  (3.20)  is 
used  as  our  closure  equation.  The  above  procedure  is  then  repeated  (if  necessary)  for  the 
final  profile  shape  using  (3.21)  as  the  closure  equation.  If  one  and  only  one  solution  leads 
to  no  contradiction,  it  is  assumed  to  be  the  true  solution.  Note  that  similar  to  closure 
method  one  where  the  possibility  of  (‘w>d'v)min  occurring  at  cloud  top  was  not  considered, 
the  case  where  ( w'9'v)h  <  0  is  not  considered.  This  is  neglected  for  the  same  reasons,  it 
adds  unnecessary  complications  to  the  program  without  providing  any  physically  realistic 
solutions. 

With  our  closure  complete,  the  same  four  step  iterative  process  described  in  section 
3.4.1  can  be  followed.  The  only  difference  is  that  (3.18),  (3.20)  and  (3.21)  now  replace 
(3.60)  and  (3.61). 

This  model  has  also  been  program  in  FORTRAN  and  will  run  interactively  on  any 
IBM-compatible  computer.  The  program,  STRATUS  Version  1.1,  requires  the  same  input 
and  produces  the  same  output  as  version  1.0. 


Chapter  4 


EXPERIMENTS  AND  RESULTS 

In  this  chapter  the  results  of  numerous  experiments  accomplished  using  STRATUS 
Versions  1.0,  1.1,  2.0,  and  2.1  are  presented.  Section  4.1  discusses  the  experimental  results 
from  the  steady-state  models  (STRATUS  Versions  1.0  and  1.1),  while  section  4.2  discusses 
the  experimental  results  from  the  time-dependent  models  (STRATUS  Versions  2.0  and 
2.1). 

4.1  STEADY-STATE  MODEL  EXPERIMENTS  AND  RESULTS 

Experimental  results  for  the  steady-state  case  are  similar  to  the  results  of  Schubert  et 
a l.  (1979)  obtained  for  similar  input  data.  The  purpose  of  this  section  is  to  examine  the 
experimental  results  over  a  wider  range  in  the  shortwave  radiation  partition,  y! .  Schubert 
et  aJ.  (1979)  allowed  y!  to  vary  in  the  range  0.0  <  y!  <  0.2.  Here,  y!  is  allowed  to  take  on 
values  as  large  as  1.0.  The  effects  of  this  change  are  examined  for  both  closure  methods. 
Also  examined  are  the  effects  of  varying  the  divergence  and  sea-surface  temperature  on 
cloud  top  and  base  heights  for  both  closure  methods. 

To  examine  the  effects  of  varying  the  shortwave  partition,  cloud  top,  cloud  base,  and 
all  convective  fluxes  have  all  been  plotted  as  a  function  of  y!  where  0  <  yl  <  1,  while  hold¬ 
ing  the  following  parameters  constant  at  the  specified  values:  sea-surface  temperature  (16 
°C),  large  scale  horizontal  divergence  (5X 10-6  s-1),  longwave  radiation  partition  (0.2),  the 
entrainment  parameter  (0.2),  and  the  location  (San  Diego,  CA)  for  both  closure  methods 
(Figs.  4.1-4.4).  Figures  4.1  and  4.2  show  the  cloud  top  and  base  height  curves  for  both 
closure  methods.  The  results  of  these  curves  are  as  one  might  expect.  As  the  shortwave 
partition  is  increased,  more  radiative  cooling  occurs  at  cloud  top  and,  as  a  result,  the 
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boundary  layer  deepens.  In  addition,  the  mixed  layer  warms  causing  a  rise  in  the  cloud 
base  height.  Both  Figs.  4.1  and  4.2  clearly  show  this  trend.  The  main  difference  between 
closure  methods  occurs  when  the  closure  equations  change  to  adjust  for  changes  in  the 
w'9'v  profile.  By  examining  the  convective  flux  curves  in  Fig.  4.3  it  is  seen  that  (w%) min 
changes  from  just  below  cloud  base  to  the  surface  at  approximately  fi'  =  0.6.  This  change 
is  accompanied  by  a  kink  in  several  of  the  fluxes.  It  is  also  at  this  point  in  Fig.  4.1  that 
a  kink  appears  in  both  the  cloud  top  and  cloud  base  curves.  This  kink  is  attributed  to 
the  change  in  closure  equations  due  to  the  change  in  the  location  of  (u7#'jmjn.  Closure 
method  two  provides  much  smoother  transitions  between  closure  equations.  In  Fig.  4.4 
it  is  seen  that  the  shape  of  the  xi/6 '  profile  changes  from  that  of  Fig.  3.2a  to  that  of 
Fig.  3.2b  at  approximately  /*'  =  0.3.  There  is  also  a  change  in  profile  shape  from  that  of 
Fig.  3.2b  to  that  of  Fig.  3.2c  at  approximately  fi'  =  1.0.  The  difference,  however,  is  that 
no  kinks  appear  in  either  the  cloud  height  or  convective  flux  curves  at  the  corresponding 
points.  Because  of  this  smoother  transition  between  closure  equations,  closure  method 
two  is  believed  to  provide  a  more  physically  realistic  solution  to  the  model. 

It  should  be  mentioned  that  although  both  models  produce  physically  reasonable 
trends,  the  solutions  for  cloud  top  and  base  height  appear  to  be  somewhat  too  low.  How¬ 
ever,  this  can  be  remedied  by  decreasing  the  large  scale  horizontal  divergence.  As  Schubert 
(1976)  noted,  halving  the  divergence  approximately  doubles  the  depth  of  the  boundary 
layer.  This  result  has  also  been  found  using  this  model,  as  will  be  shown  in  the  following 
section.  The  reason  a  divergence  value  of  5.0xl0-6  s_1  has  been  used  is  that  this  value  is 
typical  of  July  Southern  California  coastal  values  as  taken  from  the  climatological  analysis 
of  Neiburger  et  a J.  (1961). 

To  further  examine  the  effects  of  using  closure  method  one  versus  closure  method  two, 
cloud  top  and  cloud  base  height  have  been  plotted  as  a  function  of  sea-surface  temperature 
and  large-scale  divergence  for  both  closure  methods.  All  other  variables  were  held  fixed 
at  the  following  values:  wind  speed,  7  ms-1;  longwave  partition,  0.2;  shortwave  partition, 
0.2;  entrainment  parameter,  0.2;  and  location,  San  Diego,  CA  (Figs  4.5, 4.6).  Figures  4.5a 
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and  4.6a  show  similar  trends  in  that  the  cloud  top  height  increases  as  you  increase  the 
sea-surface  temperature  and  decrease  the  divergence.  Both  Figs.  4.5a  and  4.6a  agree 
with  Schubert’s  (1976)  finding  that  the  boundary  layer  depth  doubles  as  you  halve  the 
divergence.  Closure  method  two,  however,  tends  to  predict  higher  cloud  top  and  base 
values  than  does  closure  method  one  and  therefore  better  represents  the  real  atmosphere. 
Figures  4.5b  and  4.6b  show  similar  trends  for  cloud  base  height,  although  cloud  base 
height  varies  much  less  as  the  divergence  is  decreased.  Again,  however,  closure  method 
two  produces  higher  height  values  than  does  closure  method  one. 

4.2  TIME  DEPENDENT  MODEL  EXPERIMENTS  AND  RESULTS 

The  first  goal  in  producing  time  dependent  results  was  to  determine  which  closure 
method  produced  the  best  results.  This  was  done  by  inputting  identical  data  into  both 
STRATUS  Version  2.0  and  STRATUS  Version  2.1.  The  following  input  was  used:  sea- 
surface  temperature,  16°C;  divergence,  5xl0-6s-1;  wind  speed,  7  ms-1;  longwave  radi¬ 
ation  partition,  0.2;  shortwave  radiation  partition,  0.4;  entrainment  parameter,  0.2;  and 
location,  San  Diego,  CA.  The  results  of  this  experiment  were  similar  to  those  obtained 
by  the  steady-state  case.  Both  closure  methods  produced  the  same  trends  with  closure 
method  two  predicting  higher  values  of  both  cloud  top  and  cloud  base  (Figs.  4.7  and 
4.8).  The  most  notable  difference  was  in  the  u/(^  +  l')fj  curve  for  closure  method  one 
(Fig.  4.9).  It  can  be  seen  from  Fig.  4.9  that  closure  method  one  produced  a  kink  in  the 
w'(q'  +  l')u  curve  at  approximately  1030  hrs  LST.  This  kink  is  again  associated  with  the 
change  in  closure  equations  as  (tu'fl'Jmin  goes  from  just  below  cloud  base  to  the  surface. 
In  Fig.  4.10,  the  closure  equation  changes  as  the  t i/9'v  profile  changes  shape  from  that  of 
Fig.  3.2a  to  that  of  Fig.  3.2b  at  approximately  0800  hrs  LST  and  again  from  that  of  Fig. 
3.2b  to  that  of  Fig.  3.2c  at  1050  hrs  LST.  However,  the  fluxes  showed  no  kinks  at  the 
corresponding  times.  Similar  results  to  these  were  obtained  for  a  variety  of  input  data. 
For  this  reason,  closure  method  two  was  used  in  the  following  experiments. 

It  should  be  mentioned  that  some  of  the  results  in  Fig.  4.9  are  physically  unrealistic 


since  the  +  l ')w  curve  becomes  negative.  This  seems  to  imply  that  for  a  short  period 
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DIVERGENCE  (10'V) 

Figure  4.5:  Isopleths  of  (a)  cloud  top  height  and  (b)  cloud  base  height  in  meters  as  a 
function  of  sea-  surface  temperature  and  large-scale  divergence  for  closure  method  one. 
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DIVERGENCE  (lOV) 

Figure  4.6:  Isopleths  of  (a)  cloud  top  height  and  (b)  cloud  bate  height  in  meters  as  a 
function  of  sea-surface  temperature  and  large-scale  divergence  for  closure  method  two. 
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Figure  4.9:  Radiative  and  convective  fluxes  as  function  of  time  for  closure  method  one. 
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Figure  4.10:  Radiative  and  convective  fluxes  as  function  of  time  for  closure  method  two. 
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moisture  is  being  supplied  from  above  the  inversion  which  is  clearly  impossible  since  we 
assume  this  air  to  be  much  drier  than  the  mixed  layer  air.  The  problem  is  not  that  the  air 
above  the  inversion  is  more  moist  that  of  the  mixed  layer,  but  rather  it  is  due  to  a  negative 
entrainment  rate.  Because  the  shortwave  partition  is  relatively  small,  the  majority  of  the 
shortwave  heating  occurs  at  cloud  top.  Thus,  from  (3.44),  Fu  -  Fh  becomes  small  and 
possibly  negative.  From  (3.50)  it  is  seen  that  a  small  or  negative  value  of  Fu  -  Fh  will 
cause  dH/dt  to  become  negative.  If  dH/dt  is  negative  and  of  great  enough  magnitude, 
the  entrainment  rate  ( dH/dt  +  DH)  can  actually  become  negative  (which  is  physically 
unreasonable).  Thus,  from  (3.8)  it  is  seen  that  if  the  entrainment  rate  is  negative  and 
the  air  above  the  inversion  is  drier  than  that  of  the  mixed  layer  (which  is  the  case  here), 
w'(q'  +  l')H  will  become  negative.  If  we  decrease  the  shortwave  partition  slightly,  closure 
method  two  will  also  display  this  same  feature.  On  the  other  hand,  if  we  increase  the 
shortwave  partition  we  can  completely  remove  this  feature. 

The  next  set  of  experiments  looks  at  the  theoretical  effects  of  the  diurnal  cycle  on 
the  growth  and  dissipation  of  the  horizontally  homogeneous  marine  stratocumulus  cloud 
layer  as  calculated  by  STRATUS  Version  2.1.  The  modeling  results  are  compared  with 
observational  results  obtained  by  Betts  (1989)  and  Davies  and  Blaskovic  (1988). 

Betts  (1989)  completed  an  observational  study  of  this  phenomena  using  sounding 
and  ceilometer  data  collected  on  San  Nicolas  Island  off  the  southern  California  Coast  for 
the  First  ISCCP  (International  Satellite  Cloud  Climatology  Project)  Regional  Experiment 
(acronym  FIRE)  (Albrecht  et  a 1.,  1988).  The  data  represented  the  two  day  period  July 
10-12,  1987,  during  which  light  winds  (V  <  5.0  ms-1)  prevailed.  Figure  4.11  is  a  graph  of 
Betts’  (1989)  analysis  for  the  two  day  period  which  clearly  shows  a  thinning  of  the  cloud 
layer  during  the  afternoon  hours;  this  is  due  primarily  to  the  rapid  rise  of  cloud  base  depth 
and  relatively  invariant  cloud  top  height.  Similar  results  were  also  obtained  by  Davies 
and  Blaskovic  (1988)  in  which  the  diurnal  variation  in  cloud  thickness  was  observed  over 
the  whole  FIRE  experiment  (July  1-19,  1987)  using  both  columnar  liquid  water  content 
measurements  and  ceilometer  data  combined  with  sounding  data  (Figs.  4.12a,b).  These 
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Figure  4.11:  Observational  results  of  cloud  top  and  base  height  as  a  function  time  (PST). 
After  Betts  (1989). 


a.  a.  4.  a.  a.  ta.  ta.  i*.  i*.  ta.  at.  aa.  a*. 

PaaiMa  Oayli#*  Tta*  <hrl 

Figure  4.12:  Observational  values  of  (a)  column  liquid  water  as  determined  by  the 
NOAA/WPL  radiometer  and  (b)  cloud  thickness  a a  determined  by  ceilometer  and  ra¬ 
diosonde  data  as  a  function  of  time  (PST).  From  Davies  and  Blaskovic  (1988). 
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figures  indicate  that  the  cloud  thickness  reaches  a  maximum  value  of  400m  during  the 
morning  hours  and  thins  to  approximately  28m  during  the  late  afternoon. 

This  afternoon  “burn  off”  of  the  cloud  is  not  detectable  in  the  modeling  results.  It 
can  be  seen  from  Fig.  4.7  that  although  cloud  base  rises  during  the  maximum  solar 
heating  period,  the  cloud  thickness  remains  relatively  invariant  with  time  due  to  the 
diurnal  change  in  cloud  top  height.  Thus,  it  would  seem  that  the  afternoon  burn  off  of 
the  stratocumulus  clouds  is  not  valid  under  these  conditions.  However,  by  allowing  the 
majority  of  the  shortwave  radiation  to  penetrate  into  the  mixed  layer,  the  variance  of  the 
longwave  cooling  with  time  is  decreased  and  therefore  the  variance  of  cloud  top  height 
with  time  is  decreased  as  is  shown  by  equation  (3.50).  At  the  same  time,  the  variance 
with  time  of  the  mixed  layer  heating  is  increased  and  therefore  the  variance  of  cloud  base 
height  with  time  is  increased.  In  short,  we  can  maintain  a  constant  cloud  top  height 
and  allow  the  cloud  base  to  rise  simply  by  increasing  /x'.  The  end  result  is  an  afternoon 
thinning  of  the  cloud  layer.  To  demonstrate  this,  Fig.  4.7  has  been  recalculated  using  a 
fj,'  of  0.8  instead  of  0.4.  The  result  is  shown  in  Fig.  4.13. 

Figure  4.13  clearly  shows  an  afternoon  thinning  of  the  cloud  layer.  The  cloud  thick¬ 
ness  reaches  its  minimum  value  at  approximately  1400  hrs  LST  which  roughly  corresponds 
to  Betts  (1989)  and  Davies  and  Blaskovic  (1988)  observational  values.  Similar  patterns 
of  diurnal  thinning  (although  of  differing  magnitude)  are  reproduced  by  the  mixed  layer 
model  of  Hanson  and  Gruber  (1982)  in  which  the  shortwave  absorption  coefficient  is  al¬ 
lowed  to  change  based  on  the  cloud  liquid  water  content. 

Although  our  model  accurately  predicts  the  thinning  trend,  the  magnitude  of  the 
change  in  thickness  is  largely  underpredicted.  Even  if  \i!  is  increased  to  1.0  the  change 
in  thickness  still  goes  largely  underpredicted.  In  an  attempt  to  solve  this  problem,  the 
shortwave  absorption  coefficient  was  increased  from  its  value  of  7%  to  a  value  of  12.4% 
corresponding  to  a  mean  of  the  values  calculated  by  Duda  (1989)  using  data  from  July 
8  and  10,  1987,  for  the  FIRE  project.  The  specific  values  Duda  (1989)  calculated  were 
14.4±3.5  for  July  8  and  10.4±4.4  for  July  10.  The  reason  for  increasing  the  shortwave 
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radiation  absorption  coefficient  was  to  increase  the  magnitude  of  the  variability  in  the  cloud 
base  height.  The  result  of  this,  however,  was  model  output  with  no  physical  meaning.  In 
all  cases  but  one  where  a  shortwave  absorption  coefficient  of  12.4%  was  used,  w'{c{  +  t')H 
became  negative  during  the  maximum  heating  period.  The  reason  for  this  stems  from  the 
enormous  magnitude  of  the  shortwave  absorption  which  again  causes  F\j  —  Fff  to  become 
negative  during  the  maximum  solar  heating  period  for  all  cases  except  when  ji'  =  1.0. 
The  end  result  is  a  negative  entrainment  rate  which  forces  +  l')H  to  be  negative. 
In  a  further  attempt  to  increase  the  thinning  of  cloud  base,  a  shortwave  absorption  value 
of  10%  was  used.  Although  the  problem  of  a  negative  u /(q1  +  F)H  was  eliminated  for  a 
wider  range  of  input  data,  the  results  differed  little  from  those  obtained  for  a  shortwave 
absorption  value  of  7%.  This  seems  to  be  a  limitation  of  mixed  layer  models  in  general. 
In  a  higher  order  modeling  study  by  Bougeault  (1985)  which  discretizes  the  boundary 
layer  into  several  levels,  it  was  found  that  during  the  period  of  maximum  solar  heating  the 
mixed  layer  actually  became  separated  by  a  slightly  stable  region.  Observational  studies 
by  Betts  (1989)  have  also  shown  the  diurnal  separation  of  the  mixed  layer.  This  formation 
of  a  stable  layer  acts  to  turn  off  the  vapor  flux  which  decreases  the  liquid  water  content 
of  the  cloud  layer  and  thus  aids  the  diurnal  thinning  of  the  cloud.  A  mixed  layer  model, 
however,  is  unable  to  handle  such  a  separation  of  layers  and  therefore  the  thinning  of 
the  cloud  layer  is  limited  only  to  radiative  effects.  Because  of  the  problems  encountered 
when  the  shortwave  absorption  coefficient  is  increased,  the  diurnal  cycle  in  our  remaining 
experiments  is  represented  by  using  a  shortwave  absorption  coefficient  of  7%  and  a  //  of 
0.8 

One  final  feature  of  the  diurnal  cycle  that  has  been  simulated  is  its  influence  on  thin 
clouds  versus  thick  clouds.  The  case  of  a  relatively  thin  cloud  is  first  considered.  To 
produce  a  thin  cloud  the  wind  speed  of  the  previous  case  is  decreased  to  5  ms-1.  This 
decreases  the  moisture  flux  upward  from  the  sea-surface.  The  result,  shown  in  Fig.  4.14, 
is  a  cloud  layer  which  varies  from  180  m  in  thickness  during  the  morning  to  70  m  at  1500 
hrs  LST.  This  equates  to  a  diurnal  thinning  of  110  m.  These  results  are  quite  different  for 
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the  thick  cloud  case.  To  generate  a  thick  cloud,  the  sea-surface  temperature  is  increased 
to  17  °C,  the  large-scale  divergence  is  decreased  to  4.5xl0-6  s-1  and  the  wind  speed  is 
increased  to  10  ms-1.  This  has  the  effect  of  raising  the  boundary  layer  depth  by  warming 
the  mixed  layer  and  decreasing  the  entrainment  of  dry  air  from  above  the  inversion,  and 
lowering  the  cloud  base  height  by  increasing  the  flux  of  moisture  from  the  surface.  These 
results  are  shown  in  Fig.  4.15.  As  one  would  expect,  the  effect  of  the  diurnal  cycle  is  much 
less  for  the  thick  cloud.  The  thick  cloud  varies  in  thickness  by  only  less  than  50  m.  The 
reason  the  thick  cloud  shows  a  smaller  variation  in  thickness  is  that  the  flux  convergence 
of  heat  in  the  mixed  layer  is  large  due  to  the  increased  surface  temperature  so  that  the 
mixed  layer  radiative  flux  divergence  has  less  of  an  influence  on  the  change  in  9e.  The 
result  is  a  cloud  base  which  varies  little  in  time,  which  can  be  seen  from  (3.40)  and  (3.48). 

It  should  be  mentioned  that  in  several  of  the  above  results  A9e  became  negative.  This 
is  a  clear  violation  of  Lilly’s  (1968)  cloud  top  stability  condition.  However,  Randall  (1980) 
gives  a  more  accurate  stability  condition  that  includes  virtual  temperature  effects.  Under 
this  stability  condition,  the  cloud  layer  remains  stable  until  A9e  becomes  appreciably 
negative.  Typical  values  are  on  the  order  of  -2  degrees.  In  none  of  the  above  cases  did 
A9e  fall  below  this  value.  They  were  therefore  accepted  as  real  solutions. 
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Figure  4.15:  Cloud  top,  base  and  thickness  curves  as  a  function  of  time  for  tiie  thick  cloud 


Chapter  5 

SUMMARY  AND  CONCLUSIONS 

In  this  report  Lilly’s  (1968)  mixed  layer  model  has  been  generalized  to  include  the 
effects  of  penetrating  both  longwave  and  shortwave  radiation  into  the  mixed  layer  in  an 
attempt  to  simulate  the  diurnal  cycle.  In  addition  to  this,  two  closure  assumptions  have 
been  used:  one  which  is  a  weighted  mean  of  Lilly’s  maximum  and  minimum  entrainment 
hypotheses  and  one  which  considers  the  ratio  of  the  integral  buoyant  energy  dissipation 
over  the  integral  buoyant  energy  production  to  be  a  constant.  The  conservative  thermo¬ 
dynamic  variable,  equivalent  potential  energy,  was  used  instead  of  moist  static  energy  in 
order  to  help  broaden  the  audience  of  the  report. 

Modeling  studies  for  both  the  steady-state  and  time-dependent  solutions  for  both 
closure  assumptions  were  completed.  The  steady-state  results  indicated  that  both  closure 
assumptions  provided  similaf  results  for  the  same  input  data.  Closure  method  two  pro¬ 
vided  a  much  smoother  transition  between  closure  equations  which  vary  with  the  virtual 
potential  temperature  flux  profile.  Because  of  its  smoother  transition  between  closure 
equations,  it  is  believed  that  closure  method  two  provides  more  realistic  solutions  to  the 
real  environment.  It  was  also  noted  that  closure  method  two  produced  cloud  depths  that 
were  slightly  closer  to  that  of  the  real  atmosphere. 

The  time-dependent  results  were  used  to  both  intercompare  both  closure  assumptions 
as  well  as  to  investigate  the  diurnal  thinning  of  the  marine  stratocumulus  clouds.  Results 
indicated  again  that  for  the  same  entrainment  parameter,  both  closures  produce  similar 
results.  As  with  the  steady  state  case,  however,  closure  method  one  generates  kinks  in 
the  output  when  the  model  switches  from  one  virtual  potential  temperature  flux  profile 
to  another. 
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Results  from  experiments  of  the  simulation  of  the  diurnal  cycle  were  successful.  By 
allowing  the  majority  of  the  shortwave  radiation  to  penetrate  into  the  mixed  layer,  a 
diurnal  thinning  trend  similar  to  the  observational  results  of  Betts  (1988)  and  Davies 
and  Blaskovic  (1988)  could  be  produced.  Although  the  modeled  thinning  occurred  at 
the  observed  times,  it  was  much  less  in  magnitude.  In  an  attempt  to  correct  this,  a 
stronger  shortwave  absorption  coefficient  was  used.  This,  however,  resulted  in  physically 
unreasonable  results. 

The  effects  of  the  diurnal  cycle  were  also  modeled  for  relatively  thick  clouds  versus 
relatively  thin  clouds.  The  model  results  indicated  that  thin  clouds  were  much  more 
susceptible  to  diurnal  thinning  than  were  the  thick  clouds.  In  this  study  the  thin  clouds 
showed  a  decrease  in  thickness  which  was  twice  as  great  as  that  of  the  thick  clouds  for  the 
same  radiation  partition. 

Although  the  results  show  the  relative  success  of  our  mixed  layer  models,  they  also 
show  their  limitations.  The  greatest  success  comes  from  the  models’  ability  to  accurately 
simulate  the  trends  of  the  diurnal  cycle,  whereas  its  primary  limitation  lies  with  its  inability 
to  represent  a  separation  of  the  mixed  layer.  This  separation  is  caused  by  the  intense 
warming  and  consequent  stabilizing  of  the  mixed  layer  which  occurs  during  the  intense 
afternoon  solar  heating.  When  the  separation  takes  place,  the  flux  of  moisture  from  the 
ocean  surface  to  the  cloud  layer  is  effectively  shut  off.  The  end  result  is  a  thinning  of  the 

cloud  layer  which  is  much  greater  than  if  only  radiative  effects  were  considered.  It  may  be 

/ 

that  accurate  simulations  of  the  magnitude  of  the  diurnal  thinning  are  limited  to  higher 
order  models  such  as  that  of  Bougeault  ( 1985).  These  higher  order  models  contain  multiple 
levels  and  thus  can  simulate  the  diurnal  separation  of  the  mixed  layer.  Mixed  layer  models 
will,  however,  remain  an  important  academic  and  research  tool  because  of  their  simplicity 
and  computational  efficiency.  Although  their  use  may  be  somewhat  limited,  they  will  still 
be  useful  in  quickly  assessing  the  importance  of  the  predictive  variables  as  well  as  their 
influence  upon  each  other.  They  will  also  remain  useful  in  determining  the  accuracy  of 
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various  new  methods  of  paramaterizing  the  complex  physical  processes  which  govern  the 
marine  boundary  layer. 


REFERENCES 


Albrecht,  B.A.  and  R.S.  Penc,  1985:  An  observational  study  of  cloud-topped  mixed  layers. 
J.  Atmos.  Sci.,  42,  800-822. 

- ,  D.A.  Randall  and  S.  Nicholls,  1988:  Observations  of  marine  stratocumulus  during 

FIRE.  Bull.  Amer.  Meteor.  Soc.,  69,  618-626. 

Betts,  A.K.,  1989:  Diurnal  variation  of  California  coastal  stratocumulus  from  two  days  of 
boundary  layer  soundings.  Tellus,  41A,  in  press. 

Bougealt,  P.,  1985:  The  diurnal  cycle  of  the  marine  stratocumulus  layer:  A  higher  order 
model  study.  J.  Atmos.  Sci.,  42,  2826-2843. 

Cox,  S.K.,  1973:  Infrared  heating  calculations  with  a  water  vapor  pressure  broadened 
continum.  Quart.  J.  R.  Meteor.  Soc.,  99,  669-679. 

Davies,  R.  and  M.  Blaskovic,  1988:  Diurnal  variation  of  marine  stratocumulus  over  San 
Nicholas  Island  during  the  FIRE  IFO.  FIRE  Science  Team  Workshop  report,  Vail, 
CO.  July  11-15,  1988,  209-213. 

Duda,  D.P.,  1989:  Microphysical  and  radiative  properties  of  marine  stratocumulus  from 
tethered  balloon  measurements.  M.S.  Thesis,  Colorado  State  University,  Ft.  Collins, 
CO  80523.  120pp. 

Hanson,  H.P.  and  P.L.  Gruber,  1982:  Effect  of  marine  stratocumulus  clouds  on  the  ocean- 
surface  heat  budget.  J.  Atmos.  Sci.,  39,  897-908. 

Kraus,  H.  and  E.  Schaller,  1978a:  Steady-state  characteristics  aversions  capping  a  well 
mixed  boundary  layer.  Boundary-Layer  Meteor.,  14,  83-104. 

- and - ,  1978b:  A  note  on  the  closure  in  Lilly-type  inversion  models.  Tellus ,  30, 


284-288. 


60 


Lilly,  D.  K.,  1968:  Models  of  cloud-topped  mixed  layers  under  a  strong  inversion.  Quart. 
J.  R.  Meteor.  Soc.,  94,  292-309. 

Lowe,  P.R.,  1977:  An  approximating  poynomial  for  the  computation  of  saturation  vapor 
pressure.  J.  Appl.  Meteor.,  16,  100-103. 

Manabe,  S.  and  R.T.  Wetherald,  1967:  Thermal  equilibrium  of  the  atmosphere  with  a 
given  distribution  of  relative  humidity.  J.  Atmos.  Sci.,  24,  241-259. 

Neiburger,  M.,  1949:  Reflection,  absorption  and  transmission  of  isolation  by  stratus  cloud. 
J.  Met.,  0,  98-104. 

- ,  D.S.  Johnson  and  C.  Chein,  1961:  Studies  of  the  structure  of  the  atmosphere  over 

the  eastern  pacific  ocean  in  summer.  Univ.  Calif.  Pub.  Meteor.,  1,  1-94. 

Randall,  D.R.,  1980:  Conditional  Instability  of  the  First  Kind  Upside-Down.  J.  Atmos. 
Sci.,  37,  125-130. 

- ,  J.A.  Coakly,  Jr.,  C.W.  Fairall,  R.A.  Kropfli  and  D.H.  Lensschow,  1984:  Outlook 

for  research  on  subtropical  marine  stratiform  clouds.  Bull.  A mer.  Meteor.  Soc.,  65, 
1290-1301. 

- ,  J.A.  Abeles  and  T.G.  Corsetti,  1985a:  Seasonal  simulations  of  the  planetary  bound¬ 
ary  layer  and  boundary  layer  stratocumulus  clouds  with  a  general  circulation  model. 
J.  Atmos,  Sci.,  42,  641-675. 

- ,  1985b:  Key  problems  of  paramaterization  for  simple  cloud-topped  boundary  layer 

models.  A  paper  presented  at  the  Workshop  on  the  Modelling  of  the  Cloud- Topped 
Boundary  Layer  of  the  World  Climate  Research  Program,  Fort  Collins,  CO,  April, 
1985. 

Schubert,  W.  H.,  1976:  Experiments  with  Lilly’s  cloud-topped  mixed  layer  model.  J. 
Atmos.  Sci.,  33,  436-445. 

- ,  J.S.  Wakefield,  E.J.  Steiner  and  S.K.  Cox,  1979:  Marine  stratocumulus  convection. 

Part  I:  Governing  equations  and  horizontally  homogeneous  solutions.  J.  Atmos.  Sci., 
30, 1286-1307. 


61 


Slingo,  A.,  S.  Nicholls  and  J.  Schmetz,  1982:  Aircraft  Observations  of  marine  stratocu- 
mulus  during  JASIN.  Quart.  J.  R.  Meteor.  Soc.,  108,  833-856. 

U.  S.  Dept,  of  Commerce,  Environmental  Science  Services  Administration,  Environmental 
Data  Service,  1976-1980:  Climatological  Data,  National  Summary.  Govt.  Printing 
Office,  Washington,  D.C. 


APPENDIX 


In  this  appendix  the  linear  condensation  equation  (2.1)  and  its  coefficients  are  derived. 
We  begin  by  first  defining  the  equation  for  the  saturation  vapor  mixing  ratio  (?sat)  which 
is  a  function  of  both  temperature  and  pressure  as 


QiitiT  i  p)  — 


0.622e,(r) 

P  ~  e.(T)  ’ 


(A.  1) 


where  T  is  temperature,  p  is  pressure,  and  e,(T)  is  the  saturation  vapor  pressure  for  a 
given  temperature.  Since  depends  only  on  T  and  p,  it  can  be  expressed  linearly  in 
differential  form  as 

<A-2> 

Equation  (A.2)  is  a  linear  differential  equation  for  ft**  in  terms  of  temperature  and  pres¬ 
sure.  However,  since  (2.1)  is  in  terms  of  potential  temperature  9  and  physical  height  z, 
(A.2)  needs  to  also  be  expressed  in  terms  of  9  and  z.  The  conversions  of  the  differential 
of  temperature  to  the  differential  of  potential  temperature  and  the  differential  of  pressure 
to  the  differential  of  height  can  easily  be  accomplished  through  use  of  Poisson’s  equation 
and  the  hydrostatic  approximation. 

By  taking  the  logarithmic  derivative  of  Poisson’s  equation  it  can  be  shown  that 
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dT=jd9+—-dp,  (A.  3) 

9  cvp 

where  R  is  the  gas  constant  for  dry  air  and  cp  is  the  specific  heat  at  constant  pressure  for 
dry  air,  and  from  the  hydrostatic  equation,  we  can  see  directly  that 


dp  =  -pgdz, 


(AA) 


where  p  is  the  air  density  and  g  is  gravity. 


63 


By  substituting  (A.3)  and  (A.4)  into  (A.2)  while  making  use  of  the  ideal  gas  law, 
(A.2)  cam  be  written  as 

The  partial  derivatives  in  (A. 5)  are  determined  through  use  of  (A.l).  The  term 
( dq,at/dT)p  is  determined  by  differentiating  (A.l)  with  respect  to  temperature  while  hold¬ 


ing  pressure  constant.  The  result  of  doing  this  is 
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However,  from  the  Clausius-Clapeyron  equation,  we  see  that 

5  In  e,  L 
dT  ~  RVT2' 
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where  L  is  the  latent  heat  of  evaporation  amd  Rv  is  the  gas  constant  for  water  vapor. 


Substituting  (A. 7)  into  (A.6)  gives 


(dq*at\  _  Lq,at  (  P  \ 
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To  determine  the  term  ( dq,at/dp)r ,  (A.l)  is  differentiated  with  respect  to  pressure  while 


holding  temperature  constant.  The  result  is 
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Finally,  by  substituting  (A.8)  and  (A.9)  into  (A. 5)  the  desired  expressions  for  the  coeffi¬ 
cients  a  and  b  are  obtained.  These  are 


-  Lq*at  (  p  \ 

RvT9  \p-ej' 

and 

.  _  gpq,at  (  L  _  \ 
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where  the  units  for  a  and  b  are  deg-1  and  m_l  respectively. 

The  calculations  of  a  and  b  are  straightforward  with  the  exception  of  determining  e,. 

To  determine  e,  for  model  computations,  am  algorithm  written  by  Lowe  (1977)  is  used. 


